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The Horton-Strahler (HS) index r = max (i, j) + 6ij has been shown to be relevant to a number 
of physical (such at diffusion limited aggregation) geological (river networks), biological (pulmonary 
arteries, blood vessels, various species of trees) and computational (use of registers) applications. 
Here we revisit the enumeration problem of the HS index on the rooted, unlabeled, plane binary 
set of trees, and enumerate the same index on the ambilateral set of rooted, plane binary set of 
trees of n leaves. The ambilateral set is a set of trees whose elements cannot be obtained from each 
other via an arbitrary number of reflections with respect to vertical axes passing through any of 
the nodes on the tree. For the unlabeled set we give an alternate derivation to the existing exact 
solution. Extending this technique for the ambilateral set, which is described by an infinite series of 
non-linear functional equations, we are able to give a double-exponentially converging approximant 
to the generating functions in a neighborhood of their convergence circle, and derive an explicit 
asymptotic form for the number of such trees. 



I. INTRODUCTION 

Trees are ubiqutous structures which appear naturally in a large number of physical, chemical, biological and social 
phenomena, such as river networks, diffusion limited aggregation, pulmonary arteries, blood vessels and tree species, 
social organizations, decision structures, etc. They also play an important role in computer science (use of registers 
and computer languages), in graph theory, and in various methods of statistical physics such as cluster expansions 
and renormalization group. 

In spite of their apparent structural simplicity, and the large body of scientific work on trees (a sample of which 
is found in |lq| , ppf , [23| - [{27f and references therein ), they still offer challenges even related to the quantitative 
description of their topological structure. At the dawn of the science of complex networks it is therefore rather 
important to have a complete understanding of all the tree structures and their properties. 

A tree is defined as a set of points (vertices, nodes) connected with line segments (branches, or edges) such that 
there are no cycles or loops (a connected graph without cycles). For the simplest (unlabeled) rooted plane binary tree, 
each vertex has exactly three connecting branches, except for one vertex which is distinguished from all the others by 
having only two connecting branches coined as the root (R) of the tree, and a certain number of vertices with a single 
connecting branch called the 'leaves'. The height of the tree is defined by the maximum number of levels starting 
from the root (which has height 0), and it can be calculated as the maximum number of branches one has to pass 
to reach the root from its vertices (since the leaves have only one branch, it means that this longest excursion must 
start from one of the leaves). The paths from the leaves to the root define a natural direction on the tree (similarly 
to the river flow) which is always towards the levels of lower height. A tree of height k we call complete, if it has 2 k 
leaves each being a distance k from the root. 

Let us now mention three applications of the mathematics of trees which are directly connected to the so-called 
Horton-Strahler index of the tree, which is the subject of interest of the present paper. 

Originally, the Horton-Strahler index of a binary tree was introduced in the studies of natural river networks by 
Horton || and later refined by Strahler [ [l0| |, as a way of indexing real river topologies, since river networks are 
topologically similar to binary trees. By definition, a leaf has a rank of (some authors associate the value of 1), and 
a vertex has a rank of r = r(i,j) where r(i,j) is the index function with i and j being the ranks of the two connecting 
vertices from the level above. When 

r(i,j) = max(i,j) + 5 id , (1) 

the index is called the Horton-Strahler index (HS). The quantity of particular interest is the HS index of the root which 
thus categorizes the topological complexity of the whole tree. Several other quantities can be introduced in relation to 
the HS index. A segment of order k |TT| , or a stream of order k ]l2| is a maximal path of branches connecting vertices 
of HS index k, ending in a vertex with index k + 1. Let Sk{n, T) denote the number of segments of order A; of a tree T 
with n leaves, and < Lk(n, T) > is the average physical length of a segment of order k (the average < . > is taken on 
the tree T). The bifurcation ratios Bk(n,T) are defined as Bk(n,T) = Sk(n,T)/Sk+i(n,T), and the length ratios via 
Ck(n,T) — < Lk+i(n,T) > / < Lk(n,T) >. Horton and Strahler have empirically observed that for river networks 
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both the <Sfc(n) and < Lk{ri) > tend to approximate a geometric series, Bk{n) « S with 3 < $ < 5 and Ck{n) f» £ 
with 1.5 < £ < 3. Such networks are called topologically self-similar 13 1. The notion of HS index is further refined 
by introducing the biorder (i,j) of a vertex, representing the HS indices of its two children pi, O], and then 
studying the ramification matrix, with elements related to the number of vertices with a given biorder. 

Another interesting application of the mathematics of binary trees and the HS index, is in the description of the 
branched structure of diffusion-limited aggregates see Ref. (ll| and references therein. In this case the structures are 
grown on a substrate (which can be a point or a plane) by letting small particles diffuse towards the aggregate where 
they stick indefinitely at their point of first contact with the cluster. This creates complex and involved branched 
structures, whose topological complexity still remains a challenging problem to describe. 

Finally, the last application we would like to mention is known as the word bracketing problem j^] which has 
obvious implications in computer science. Let us consider an alphabet of n letters, A = {Y±, Y2, Y n } and a word 
S = XiX2X3..x n —iX n , Xi S A. A 2-bracketing of the word S is a partition of its letters (by keeping their order) in 
groups of two units enclosed in brackets, where a unit can be a letter or a subpartition enclosed in brackets, such 
as (xiX2){x3(x4 : X5)) , or (xi(x2(x3(x4,X5))), etc. The bracket (1*1^2) between two units may be associated with a 
multiplicative composition law (U1U2) = u\ ■ U2 G A. For example let the alphabet A be all the positive integers, 
and the composition law be the regular multiplication of numbers. Then a bracketing of the multiple product S 
corresponds to one particular way of calculating S. A one-to-one correspondence can be made immediately to trees: 
let the letters x\, X2, x n of the word S be associated with the leaves of a binary tree. To a particular bracketing 
of S it corresponds a particular tree constructed by associating a lower level vertex to a bracketing (111U2) (one may 
think of the brackets as representing the branches of the tree). The main question is how many ways are there to 
calculate such a product. If one assumes that the multiplication law is neither associative nor commutative, then the 
problem is refered to as the Catalan problem, see Ref. || for a number of solutions. The number of such bracketings 



is given by the Catalan numbers, a r , 
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. The corresponding set of trees (see Fig.l for n = 4) is in fact the set 



of rooted, unlabeled binary plane trees according to this bijection. 




2) 





5) 




FIG. 1. The set of rooted, unlabeled binary plane trees corresponding to all the possible non-commutative, non-associative 
bracketings of the four letter word abed, n — 4. 

For later reference, we mention that the generating function A(£) — X)^°=o£™ a ™ °f the Catalan numbers obey the 
equation A 2 - A + £ = 0, with A(0) = 0, so 



A(S) = -{l-y/l=4Q. (2) 

The power series A(£) converges within a disk of radius a c = 1/4. 

The problem of enumerating trees becomes more difficult if the composition law is commutative, which was first 
studied by Wedderburn and Etherington (WE) ||, jj, In the tree language, this means that two trees are 
considered identical if after a number of successive reflections with respect to the vertical axes passing through the 
vertices they can be transformed into each other and in this case they are said to be homeomorphic [Q . For the example 
shown in Fig. 1, there are only two such trees, since trees 1), 2), 4) and 5) can be transformed into each other. The 
trees that cannot be transformed into each other are called non-homeomorphic. The set of non-homeomorphic trees 
is called the set of ambilateral trees, jll|, p% . Let the number of such trees with n leaves be denoted by w n . The 
generating function (GF) defined as IF(£) = X^^Lo^™ 10 ™ obeys the nonlinear functional equation: 

w(o = z+^w(o 2 + \w(e) (3) 

which has extensively been studied by Wedderburn H. Otter B studying a more general counting problem where 
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the vertices can have at most m branches comes to the conclusion that for the ambilateral trees, if n is large we have: 



c n 



-3/2 n 



7" (4) 



where 7 = 2.4832535.... The method developed by Otter gives an iterative approach to c and 7. For example 7 is 

7 = lim 4 " (5) 

n — >oc 

where sq — 2, s n — 2 + s^ l _ 1 so that for n = 4 one already obtains an extremely close value of 7 ~ (2090918) 2 . 
Later, Bender developed a more general approach Jl6| deriving the same results. The coefficient c in (Q) can also be 
computed: c = 0.31877662.... 

The more practical application of the bracketing problem within computer science is the computation of arithmetic 
expressions by a computer. A general arithmetic expression involving only binary operators can simply be mapped 
onto a binary tree, called the syntax tree, which has as leaves the operands and the inner vertices the operators. A 
computer traverses this tree from the leaves towards the root and it uses registers to store the intermediate results. 
In general there are many ways of traversing such a tree, and the program that uses the minimal number of registers 
is the most efficient, or optimal one. Ershov has shown |l7j ] that the optimal code will use exactly as many registers 
to store the intermediate results as the HS index of the associated syntax tree. 

In the present paper we investigate how the HS index is distributed on both the rooted, unlabeled, plane binary 
set of trees, and on the ambilateral set of binary trees. We first answer this question on the rooted, unlabeled, plane 
binary set, since it is simpler, but it will also provide us with a technique that can be extended to tackle the problem 
for the ambilateral set. For this set, the question was first answered by Flajolet, Raoult and Vuillemin, jl8| with a 
method somewhat similar to the one presented here. The enumeration problem of the HS index on the ambilateral set 
is, however, inherently more difficult since it involves functional equations with nonlinear dependence in the argument 
similar to Eq. (||) , and therefore an explicit solution in a closed form becomes impossible to attain. The derivation of 
an approximant formula for the number of ambilateral trees sharing the same HS index at the root is the main result 
of this paper. 

The paper is organized as follows: first we present our derivation of the enumeration problem for the HS index on 
the unlabeled set in Section II, and then use this method of derivation from this case to develop a technique that can 
be used to attack the enumeration problem on the ambilateral set in the asymptotic limit, presented in Section III. 
Section IV is devoted to conclusions and outlook. 



II. DISTRIBUTION OF THE HS INDEX ON THE UNLABELED SET 



Let us observe that the root R of the tree has always two subtrees attached to it via the two branches, with k and 

(r) 

n — k leaves, respectively, k = 1, 2, .., n. Let Nn ' denote the number of unlabeled trees with n leaves that share the 
same HS index r at the root. A recursion is found for this number in the light of the observation above: 

n— 1 I r—1 r—1 1 

K r) - E ^t 1 ^ + K r) E E < j) (6) 

fe=l [ ]=0 j=0 J 

with the conventions = 0, Nn = $n,i, N^ = S r .o- If the generating function for the variable n is defined as 

D A0 = EZo^^' then it obe y s: 

r-1 

D r = D 2 r _ 1 + 2D r y £,D j , r>l, D = £ (7) 

3=0 

Next we give an exact solution to (0). Let us introduce the sum B r = X^=o Dj> r ^ L B = 0, B\ = £ . Then 
D r = B r+ i — B r , and after rearranging the terms, Eq. (Q) becomes G r = G r _i, where G r — B^ + B r+ i{\ — 2B r ). 
This means that G r = Go = £, i.e.,: 

B 2 r + B r+1 (l - 2B r ) = £ , r>0 (8) 

Note that the left hand side of (||) remains invariant to B r i— > 1 — B r which is another solution of (j^). However, since 
in case of the HS index Bq = 0, this latter solution has to be dropped. If we make 2B r = 1 — C r , (||) simplifies to 
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C 2 — 2C r C r +i = 4£ — 1. Which, after dividing on both sides by 4£ — 1, and introducing Z r = C r /V4£ — 1, becomes: 

z -» = ^ < 9 » 

Let us now write Z r = cot(tv), such that vq — arctan (-\/4| — 1). Then @ becomes cot(iv+i) = cot(2u r ) which leads 
to = 2v r + 7rm, raeZ, and which in turn is solved easily. Thus, Z r = cot(2 r «o), so one finally obtains: 

A.(0 = - . -— . (io) 

2 sin (2 r+i arctgV4£ — lj 

Eq. (10) is the exact solution to (^) in the complex £ plane. On the real axis, within the radius of convergence a c 
the above expression takes the form: D r (f) = VI -4£/ [2sh (2 r+1 arcthVl - 4f)] , £ < a c = 1/4. Since within the 
convergence disk one must have J2^=o -^riO = we just obtained the identity (using (0)): 

1 +E^) =cth ^ * >0 c n > 

r— 1 v 7 

This identity can be checked to hold via more direct methods p2|. The singularities of D r (^) lie on the positive real 
axis at: 

finity (cc 

if one simply iterates (0) we obtain 



with an additional singularity at infinity (corresponding to k = 2 r ). We certainly have £[ r ' > a c . On the other hand 



where P r (£) is a polynomial in £ of order 2 r - 1: P^) = - £, P 2 (£) = Pi(£)( 2_1 - 2 £ + £ 2 ), P 3 (£) = P 2 (£)( 2 ~ 1 - 
4£ + 10£ 2 — 8£ 3 + £ 4 ),... etc. One can find an explicit form for this polynomial from the general solution ( |To| ) if 
one invokes the identity sin (nx) — nsuixcosxYl^i 2 ^ 2 (l — sin 2 xj sin 2 (A;7r/n)) , so that ( |T3| ) is recovered with: 

Pr(£) = Tlk=i ct g 2 (^/ 2r+1 )(£fe — £)• ^ is eas y to show, however, that IljLi 1 ctg 2 (fc7r/2' +1 ) = 1, so the polynomial 
simplifies to: 

2 r -l 

fc=l 

expression valid on the whole complex £ plane. Based on the explicit solution we obtained, one can give an exact 
form to the distribution of the HS index on the unlabeled set of trees, by inverting the generating function via: 



^-Si/^o-sj/^ffe) (15) 



One can write: 



2 r -l j^ir) 



2 rp r{0 (16) 

where — 2~ r Yl ~ > 3 = 1, 2 r — 1. By Cauchy's theorem the integrals in (|l^) are readily 

k = l \ ' 

performed, and one obtains: 



- - j > r (r) -| -(n-2 r +l) 
^3 



3 = 1 



> n> ~ T (17) 
< n < 2 r - 1 
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From (|T|) it follows that Y?j=i I^P = 2 ~ r / P r(0) = 1. To obtain the last equality we used the form Q and 



(|12j). Thus: N$) = 1, r = 1, 2, ... The numbers AJ> can be calculated as follows: observe that 



|(r) 



AV = liui 



(18) 



where we used the L'Hopital rule in the last equality. On the other hand from ( |l3| ) and ([!(]) it follows: 2 r P r (£) = 
2£ 2r sin (2 r+1 arctg^/4£ — l) /\/4£ — 1. Taking the derivative of this equation at the point Q , and inserting it in Ji"s| ) 
it yields: 



,4 



(r) 



4£ M - 1 



(-1) 



2r+l^( r ']2 r -l 



Thus, we obtain from ( |T^ ) the more explicit form 



or using 



(19) 



1 2r - x 4f (r) - 1 

" 2 r+1 ^ [£ ]™ ~~ 



(20) 



7V (r) 



4" 



i=i 



,1 + 1 ■ 2 I 3* 

2 r + 1 



J7T 
2 r+l 



(21) 



an expression first derived by Flajolet et. al. |18[ ]. Following this paper Jig ], our polynomials P r can be simply 
connected to the Tchebycheff polynomial U Q, via the relation: P r (f) = 2-''£ 2r ~ 1 / 2 [/ 2 ,.+i-i(l/(2\/£))- 

If one employs the Poisson resummation formula for functions defined on a compact support (see Ap- 



pendix B in Ref. |l9[) on 

E^iV 2 [( r -)](fc)|^ 

method, see [fL8f . 



2l|), an equivalent combinatorial expression can be derived in the form: N 1 



(r) 



l+(2m-l)2'' 



where (V/)(fc) = /(fc) — /(fc — 1) is the finite difference operator. For a different 



(r) 

Scaling limits. Next we briefly present the results of an asymptotic analysis on the Nn numbers. Since N, 
is an enumeration result, it typically contains several scaling limits. In physical processes, during the growth of 
branched structures, usually only one of these limits is selected, and in frequent cases this limit has self similar 
properties (such as for DLA, or for random generation of binary trees, |2(|). By definition, the family of trees that 
obey Um r _ ) , 00 (In [n(r)])/r = const. = \nB is called topologically self similar jl3[ ], where B is the bifurcation number. 
1) n — > oo and r fixed. In this case the first term in (Tu) dominates the sum and the asymptotic behavior is given by 



(r) 



N, 



(r) 



a nln [4cos' i (7r/2 r+1 )] 



The rate of the exponential growth is a number between In 2 and 2 In 2. 



2 

2) n — * oo, r — > oo, -y/n/2 

corrections) and yields: Nn^ ~ 7r 2 2~ 3 ( r+1 ) e "( 21n2_7r /4 J. If v / n/2 r diverges with r slower than exponential, we 
have topological self similarity with B — 4. 

3) n — > oo, r — > oo, y/n/2 r — > g?, with some < d < oo. In this case the rest of the terms in ( [H| ) (after the first has 



oo. Here the first term in ( pT| ) is still dominant (the rest being exponentially small 
" 2 /4 r+ Tlf 



been factored out) are of the type j x 



T (r) 



and the final expression is: 



A(d) 4 n n- 3 / 2 4 n . The topological 
self similarity is obvious with B = 4. The factor A(d) is given by A(d) = ir 2 e~* 2d2 (1 - e _7r2rf2 )/(l + e _7r2d2 ) 3 . 
4)rt — * oo, r — ► oo, y/n/2 r — ► 0, and n/2 r — > oo. In this case the analysis is performed easier from the combinatorial 

(r) 

expression of N„ and yields: N, 



-1/2 „-5/2 n21n2-4 r /rs 
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III. DISTRIBUTION OF THE HS INDEX ON THE AMBILATERAL SET 



Let us now analyze the same question on the set of ambilateral trees, and denote the number of ambilateral trees 



Cr) 

with n leaves and HS index r by AfA ■ We certainly must have the relation 

oo 
r=0 



(22) 
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FIG. 2. Particular values for the number of ambilateral trees 
there is no such tree. 



with n leaves and HS index r. The shaded entries mean that 



The table in Fi. 



Mn^ = 6 n ,i, and 



ig. [2] gives the distril 

MP = 1 - 8n,U SO 



distribution of the HS index for n up to 32 and r — 2, 3, 4, 5. We can check easily that 
for simplicity these are not represented in the table. 

The numbers Al„ ' obey slightly more complicated recurrence relations since now the counting has to be done on 
a more restricted set. We must distinguish between odd and even n values. However, the two cases can be combined 



into one, if the convention M, 



(r) 



0<k<j<n 



for v non-integer is adopted. The corresponding recurrence relation becomes 

1 



r-l 

•E 

s=0 



M 



(«) 



M[ s) M^ r) 



M 



(r) X- M W 



s=0 



/2 



-M 



(r-l)i 
i/2 1 



1 



n/2 , 



(23) 



G 



The generating function K(0 = E^o Mn' > wm thus obey: 



1 [V r -l(0] 2 + Vr-l(e) 

2 i - E K(0 



VAO = ^ 1 > r > 1, (24) 



and Vo(£) = £• As a check for the correctness of (p4|), let us see if we recover the identity E^oUr(£) = 
(which follows from (||)). Eq. (||) is equivalent to 2V r {£) - 2Y? s Zl V a (£)V r (£;) = [V r -i{£)] 2 + K-i(£ 2 )- Introduce 
the temporary variable G(£) = E^-p V r ($) and sum both sides of the equation over r, r — 1,2, ..,oo. One obtains 

2(G(£)-0-2E~i£llo^)^(0 = £~ [Vr(0] 2 +G(e)- Using the identity 2 E~ i E^o ^ ^ = (E^o^) 2 - 
E^LoK 2 , one finds G(£) = £ + ±[G(£)] 2 + ^G(^ 2 ) which is precisely Eq. (|), showing that G(f) = W(Q, i-e., the 
relation E^lo = W(0 holds, indeed. 

In contrast to the previous case, the functional recurrence ( p4| ) cannot be treated in an exact analytical fashion due 
to the functional dependence on £ 2 . However, one can derive the asymptotic behaviour and make statements that 

(r) 

will lead to rather close approximations of the M„ numbers. It is instructive to look at a few particular values, first: 



e 



(1 - 20(1 - 2£2)(1 - C 4 )(l - 3£ + 4^ 3 - C 4 ) 



Inverting V^O, one obtains: = [2"- 1 - 3 + (-l) n " 4 ]/6, n > 4, which can be checked to hold, sec the table in 

Fig. ||. The result from the inversion of V3 (£) is already so complicated that it is not worth presenting. As the index 
r increases, the polynomial expressions become more and more involved. Figure ^| shows the function Vs(£) in the 
interval [-2.0,2.0]. 



10 



5 







-5 







I 


\ 


1 "us 


a 

|P 


1 

f\ 
II 




! J 






t 


i ■ p!?~ 

f 






1,1 j 


1 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 



FIG. 3. The generating function Vs(£) on the real axis. The function was evaluated in more than 1.3 • 10 5 points, and 
represented by dots. 

For every r, the power series for V r (£) has non-negative coefficients, > 0. Based on a classic theorem of 

complex analysis, this means that on the circle of convergence, of radius a r > 0, there will be a singularity of V r (£) 
at £ = a r . Next we show, that we have the ordering < a r +i < a r < 1 for r > 2, and the limit lim^oo a r exists 
and it is equal to a = 1/7 = 0.4026975036.... We shall use mathematical induction to prove the ordering. From the 
particular examples above it follows that 02 = 0.5, a.3 — 0.424507.... Let us now assume that oy < ctj-i < 1 for all 
j < j > 2- We will show that a r+ i < a r . Note that the radius of convergence for Kj(£ 2 ) is ^fa] > otj, if atj is less 



7 



than unity. By reductio ad absurdum, let us assume first, that a r +\ > a r . This means, that is analytic in a r . 

From (U), 

By the argument above, V r (£, 2 ) is analytic in a r (its radius of convergence is ^/oh > ctr, since by assumption 
a,. < a r _i < ... < a2 = 1/2 < 1). In the denominator of (p6|), aZZ functions Vj, j = 0, 1, ..r — 1 are analytic in a r , 
because by assumption they all have radii of convergence strictly larger than a r . However, V r is singular in a r , and the 
singularities do not cancel in the numerator and denominator of (|26[), and thus V r +i is singular in a r , a contradiction. 
We are left to prove that a r+ i — a r cannot hold. Let us denote B r = 2 s =o Again, we assume, that a r+ i — a r is 
true. It is easy to show, that for any finite r, |V^(av)| = oo. This means from the recurrence relation that 

B r _i(or) = l (27) 

(in the numerator of (53) we have only functions analytic at a r ). Since V r+ i(^) — + ^r(£ 2 )] /[I — Br(0]i from 

the assumption a r +i = a r it would follow that the equation B r (x) = 1 cannot have any solutions (V^+i is analytic 
within the circle of convergence) in the interval < x < a r . (Note that in the interval < x < a r , the numerator 
Vr{£) + K(C 2 ) cannot be zero, since the power series V r has only positive coefficients). The equation B r (x) = 1 is 
equivalent to B r -i(x) + V r (x) = 1. However, from 1 - B r -i(x) = [V^ix) + V r (x 2 )] /V r (x), Thus, the equation 

V r 2 (x) = V^ix) + V r -i{x 2 ) (28) 

should have no solution in < x < a r . If x is arbitrarily close to a r , then V 2 {x) is arbitrarily large. However, 
since a r ~i > a r , V r r 2 „ 1 (a:), and V r r 2 _ 1 (a; 2 ) are both bounded from above. Thus, for x sufficiently close to a r , we have 
V 2 (x) > V 2 _ 1 (x) + V r -i(x 2 ) . On the other hand, the HS index of a tree T equals to the height of the largest, complete, 

(r) a (r) 

balanced tree embedded in T. This means, that M„ = for n = 0, 1, 2, .., 2 r — 1. Also, M% r = 1. In other words, 
one must have V r (x) = x 2 (1 + 0(x)). 




FIG. 4. A magnified region of Figure |4| The arrows indicate the positions a = 0.40269... and ^fa — 0.63458... on the real axis. 

This means that V^x) = x 2 " (l + 0(x)), V r -i{x 2 ) = x 2 '' (l + 0{x 2 )), and V 2 (x) = x r x r (l + 0(x)). Since 
V 2 _ 1 (x) + V r -i(x 2 ) = 2x 2 (1 + 0(x)), there will always be an x > 0, (x < 1), sufficiently close to zero, such that 
V 2 (x) < V 2 _ 1 {x) + V r ^i{x 2 ). Therefore, there must exist an < x < a r , for which ( pq) holds, which is a contradiction. 
Thus, we have proven that < a r+ \ < a r < 1, for all r > 2. As a matter of fact we have also shown, that the radii 
of convergence satify: 

V 2 (a r+1 ) = V^iar+x) + V r ^(a 2 r+1 ), r > 1. (29) 
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Since the series a r is monotonically decreasing, and bounded from below, the limit a = lmv^oo a r exists. 

We have shown that J^^Lo ^r(£) = Since the radius of convergence for the left hand side is the minimum 

of all the radii of the terms in the summation, i.e., a, it must equal to the radius of convergence for W(£), which, as 
shown by Otter and Bender is 1/7, lim r _» 00 a r — a = 1/7 = 0.4026975036.... Taking the limit r — > 00 in (]27]), we get 

W{a) = 1 (30) 

(since by definition B r = J^^Lo Vs, so linir— >oc B r — W(£)). Or, using (||): 

W(a 2 ) = 1 - 2a (31) 

an identity also shown by Bender. Eqs. (|31) and ([?]) can simply be combined to give the iterative computation of a 
in the form already mentioned in the Introduction, as follows: if we make the temporary notation 

[7(0 = [1 - W(£)]/y/i . (32) 

Eq. (H) takes the form 

£/(£ 2 ) = 2 + [/ 2 (£) , (33) 



and Eq. (31) simplifies to 

U{a 2 ) = 2 . (34) 

Let S{x) = 2 + x 2 . Then, from © f/(£ 2 ) = S(U(0), or U(£) = S^tf 1 /*)) = S^U^*))) = S(...S([/(£ 2 -"))...), 
where there are a total of n compositions for the S function, n arbitrary. Let us now choose £ = a 2 . This means, 
[/( a 2" +1 ) = 5(...5(?7(a 2 ))...) = S(...(S{2))...), by virtue of (§|). From ©, C/(a 2 " +1 ) = [1 - T4^(a 2 " +1 )]/a 2 " . We 
have shown previously, that a < 1 (it is the limit of the monotonically decreasing series a r < 1), therefore we have: 

a = lim ( = lim s - 2 - (35 ) 



n — >oo 



since W(a 2 " +1 ) — > W(0) = 0, and where s n = S(s n -i), s — 2, just as in the Introduction. The convergence is 
double-exponential, very fast. 

(r) 

As in Section II, the asymptotic behavior of the Mn numbers for relatively large n and r is governed by the 
innermost singularity of V r (£) on the real axis. The graph of Vg shown in Figure suggests, that the generating 
function is in fact well behaved in a certain interval to the right of the radius of convergence, as, see also Figure 
^. The existence of this interval comes from the fact that the singularities of the term with nonlinear argument 
Vv_i(£ 2 ) in the numerator of (|2~I| ) kick in only beyond the circle of convergence of Vr^i^ 2 ), which is y/a r -i > a r -\. 
Thus, in the interval a r < x < ^Ja r ~\ the term with the nonlinear argument is analytic, which ultimately is 
responsible for this nice behaviour. Because, a < a r _i, for convenience we shall define the interval of this nice 
behaviour to be I — [a, y/a). In order to exploit this observation, we shall first rewrite the recurrence relation j24|). 
Let us denote H r (£) = 1 — X)I=o^«(£)- With this notation, ( p4] ) takes the form G r (£) = Gr-i(Q, r > 1, where 
G r (£) = H 2 {C) - 2H r {£ y )H r+l {£ > ) + H r (£ 2 ). This leads to the new recurrence: 

H 2 (0 - 2JT r (0flr+i(0 + H r (e) = 2£ , (36) 

Hq(£) = 1. This would be exactly solvable if it were not for the dependence on the nonlinear argument £ 2 . Note 
the resemblance to (||). Let h r (£) — 2£ — i/ r (£ 2 ), which is an analytic function in /. We also have A/i r (£) = 
h r (£) — h r -i(£) — V r (t; 2 ) = £ 2 (l + 0(£ 2 )), the latter equality being shown previously. This shows, that in the 
interval /, the r-dependence weakens extremely fast, double-exponentially with increasing r. As a matter of fact, an 
upper estimate is 

Ah r (0 < a 2 "' 1 . (37) 

In particular, Ah 3 (g) < 0.0263, A/i 4 (£) < 0.0006916, Ah B (g) < 4.79-10~ 7 , Ah 6 (0 < 2.28-10- 13 , A/i 7 (f) < 5.22-10- 26 , 
etc. Therefore, from the point of view of the asymptotic behavior, the h r functions can be replaced by their asymptotic 
expression (as r — > 00): 

HO = W(e) +2£-l. (38) 







Figure || shows the functions h r on the interval I for r = 1, 2, 3, 4, 5, 6. 

2 



1.5 

1 



0.5 







-0.5 

0.35 0.4 0.45 0.5 0.55 0.6 _0.65 

a £ \fa 

FIG. 5. The functions /i r (£) are analytic on J. This figure shows h r (£) for r = 1, 2, 3, 4, 5, 6. The convergence on I to is 
double-exponentially fast. The thick vertical lines delimit the edges of the interval I. Close to a, the h r functions cannot be 
distinguished on / for r > 3. To the right from / the h r functions develop singularities. The point y/a is a left accumulation 
point for the series of the leftmost singularities of h r (£) as r — > oo. 




Thus, instead of Eq. (|36|) we will consider: 



H 2 r (0-2H r (0H r+1 (0 = h(0. (39) 



The recurrence (p9J) in turn is easily solved in the way shown in Section I. The result is: 



H r (0 = VW) ctgl 2'-'*arctg l^jf^ J J (40) 

where ro for the moment is an arbitrary (positive integer) index. Recurrence ( |39| ) will become a good approximation 
to the recurrence ( |36| ) from an index r on. The larger r is the more accurate the approximation. Recurrence (39|) 
is applied then with initial condition iT ro (£) = H ro (£), which for modest ro values can be obtained by iterating (3q) 
ro times. 

What is the error we make when one replaces /i ro (£) with /i(£) on II Summing the differences ( |37| ) from r + 1 to 
infinity, one obtains the estimate: h(£) — h ro (C) ^ o* 2 X/m=o 

a r°(2 m -i) K -^ rw . Thus, for example, h(£) - h 5 (Z) 

is smaller than 10~ 7 , h(£) — hg(^) is smaller than 10 -13 , etc. 
Therefore, we can finally write on /: 

m} vy ?HK-i T>n> _ (41) 



In Fig. (I^) we plot the rhs of (|4l|) and the function from iterating (g4J). Note that the approximation is very 
good, and it becomes virtually indistinguishable from the true function the closer £ is to a. Larger ro values will 
also give better approximations, since the approximation is only applied from the ro index on. However, ro cannot 
be taken too high for approximation purposes, since it assumes that the exact expression of H ro (or V ra ) is known. 
This makes only the modest ro values (less than 5) useful. On the other hand, expression (^l]) is very practical in 
analysing the singularities of V and give rather close approximant expressions to these singularities. In particular, 
we see that within the interval /, ( f4l|) preserves the property that if cv is a singularity of V r i (or a zero of H r t) 
then it is a singularity of V r (or a zero of H r ), whenever r > r' . If one is interested in the asymptotic behavior, then 
a more tractable expression can be derived for the rhs of (pd]): the function /i(£) is analytic on the interval /, and 
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since already for modest r values, the innermost singularity of V r (denoted a r ) is extremely close to a, one can safely 
replace h(£) in this neighborhood by: h(£) ~ h'(a)(£ — a). 



v 




FIG. 6. The true V r (£) function (dashed line) from iterating (^), and the approximation in Eq. ( ftl] ) (solid line) for r — 8 
with ro = 3, and r = 3 with ro = 2 (the inset). 



This leads to the approximant: 

V r (0 ~ tf P (0 



sin (2 r+1 arctg (^-y/^o)) 
for sufficiently large r (here "large" means r > 4) where 



/i = \Jh'(a) , 



2 r °H ro (a) ' 

Next, we compute h f (a). One can use a very similar method to the one employed to obtain (35), to give: 

fin 



h'(a) 



lim 

n^oo S Sl...S n -l 



3.1710556. 



(42) 



(43) 



(44) 



2a 3 



so, fi = 1.780745815.... If one computes for ro = 3, we have Hj(a) = (1 — 3a + 4a 3 — a 4 )/(l — 2a — < 
0.164518.., and thus 6 = 1.3530022.... If we were to use r = 4, then one would obtain H 4 (a) = 0.082262, so 
8 = 1.3529529245 and slightly improve the approximation on 9. No significant improvement will be obtained with 
larger ro values. Figure M shows the agreement of the form given in (}42|). For clarity, we defined the function f(z) 
given by: 



tg(^) 



(45) 



Here we use the true V r function using numerical iteration of (J2J), and evaluate it in the points £ = a + 9~ 2 tg 2 (^ttt)- 
If the approximation (E^) is good, then one should have f(z) — sin(z). As seen from Fig. ^| the approximation is 
already excellent for r — 4 close to a (which corresponds to the z — point in this plot). The interval / in these 
transformed corrdinates corresponds to (0, 2 r+1 arctg(#-\/ ^fa — a)) = (0,0.577435486 • 2 r+1 ). There are no fitting 
parameters, we used for /i and 9 the values derived above. 

(r) 

In order to obtain the approximation to the number MA of ambilateral trees with the same HS index at the root, 



we will have to invert 



The singularities of the rhs of (|42j) are given by: 



kit 
T+ 1 



k = 1,2,3, ...,2 r -l 



(46) 
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(at the moment we do not care whether some of the singularities will fall outside the interval /, we just simply 
want to invert and then at the end keep only those terms from the final expression that were generated by the 
singularities within /). 




FIG. 7. The goodness of (j42|). For fi and 6 we used the values derived in the text. 
In a similar manner to the previous section, we first bring K r to an inverted polynomial form: 

1 2 r 



2 r + i r +1 - 1 Q r (o 



(47) 



where Q r is the polynomial: Q r (£) = II^—i 1 Kfc^ — ^) * The case f rom the previous Section II corresponds to 



fj, = 1, 6 = 2 and a — 1/4. Thus, if we denote by M n the numbers coming from the inversion of K r (£), then: 

- {r) _ n 1 i d£ [l + 9HZ-a)f 



We have: 



-(0 



3 =1 Cj S 

After performing the integrals, one obtains: 



r4 1 with = ; — r 

(r) _ f ' J 2 r + 1 Q 2 



E 



2 r -l 

in 



=1 Sfc 



(48) 



(49) 



min{n,2 r } 



M^Ewr-- 1 E 

j — 1 m— 



(l-ar) 



2l2 r — ro 



Q2<r( r ) 



(50) 



— ~( r ) (r) 

This expression shows that the M n may only approximate the M„ numbers in a certain limit. This is seen from 
the fact that while one must have = for n < 2 r , and M^) = 1, this is not respected by ( |50| ) (it would 

only be respected if a — 9~ 2 , however, this is not the case, and the reason behind this discrepancy is the neglected 
nonlinearity from the calculations). The limit, in which the approximation becomes good is for r large (it means 
r > 4) and n 3> 2 r . In this case the sum over m can be performed, and one obtains: 



M 



{r) - E^r^r^ 1 



a 



2 r 



(51) 
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The numbers can be calculated in exactly the same way we did in the previous section. This leads to: 



Af = (-iy +1 - 



a 



2 r -l 



(52) 



Inserting it into it yields: 



M 



yr a v I 



+» 3 («r-«)K«r-°) 



Wln+l 



(53) 



As a check to the correctness of ( |53| ) we can take = 1, = 2 and a = 1/4 from the unlabeled case, to obtain (|2C|). 
Equation ([53] ) explicitely shows the contribution of each singularity. However, if we want to approximate the Aff' 
numbers, we should also account for the condition < y/a. Using the expression (46), this leads to j < J r , where: 

2 r+l 



arctg(0v/V«- a) - (0.1838035250..) • 2 



r+l 



(54) 



Thus, using again (46): 



{r) „ W ) + i tan 2 (# £ r)[l+tan 2 (# £ r)] 

" »H i [ a + ^ 2tan 2 ( _2 ?r) ]«+ 1 



(55) 



When the asymptotic limit is generated by the innermost root a r ~ £i , i- e -i by the first term in (£>5|), one obtains 
for the topologically self similar ambilateral trees, the scaling behaviour: 



-(r) 2/i7T 2 (i 3 



Ml' ^r^e"^n- 3/ V 



*(9 3 



(56) 



and therefore B = j = l/a = 2.4832535... 

Let us now see how well formula (|5J|) approximates the M^ r ' numbers. To do this, we shall define the error 
Q { n ] = \\M ( ^ - M^ r) \/Mi r) ] ■ 100%. For example from the Table in Fig. | M 3 4 2 = 413083691. The formula above 



gives [M^] = 445781858, and thus = 7.915...%. Further error values: Q^o = 5.34132...%, Q$ Q = 0.05391...% 
= 0.003551...%. 



IV. CONCLUSIONS AND OUTLOOK 



Combinatorial enumeration of trees is typically difficult to solve when the set under enumeration obeys symmetry- 
exclusion principles, such as for the ambilateral case treated here. These symmetry-based constraints may arrisc 
in realistic situations and thus forces us to enumerate classes of subsets of trees. In the ambilateral case a class is 
defined as being formed by those binary trees that have the same number of leaves and HS index at the root and 
can be obtained one from another via successive reflections with respect to the nodes of the tree. Certainly, the 
symmetry operation defining the class must be an invariant transformation of the topological index (HS in our case) . 
An other example of such symmetry-operation-generated class-enumeration is the case of the "leftist trees" playing 
an important role in the representation of priority queues, first shown by Crane p3[ , followed by Knuth pH , who 
gives their explicit definition. An elegant enumeration for the leftist trees, using generating function formalism was 
only given very recently by Nogueira p5| . 

The existing solutions to such class-enumerations on trees (such as ours and that of Flajolet et. al. |i| and of 
Nogueira [p5[) are obtained via methods taylored for the particularities of the set and symmetry operation in question. 
It is desirable to have, however, at least on a formal level, a general encompassing theory of class-enumerations of 
topological indices. In this direction, powerful methods such as that of the antilexicographic order method developed 
by Erdos and Szekely p6| , or the method of bijection to Schroder trees developed by Chen may turn to be effective 
after a suitable extension to include topological indices such as the Horthon-Strahler index. This, however, stands as 
an open problem. 
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